Spontaneous Imbibition Dynamics of an n-Alkane in Nanopores: 
Evidence of Meniscus Freezing and Monolayer Sticking 



Simon Gruener and Patrick Huber 

Faculty of Physics and Mechatronics Engineering, 
Saarland University, D-66O4I Saarbriicken (Germany) 
(Dated: October 22, 2009) 

Capillary filling dynamics of liquid n-tetracosane (n-C24H5o) in a network of cylindrical pores 
with 7 and 10 nm mean diameter in monolithic silica glass (Vycor) exhibit an abrupt temperature- 
slope change at Tg = 54 °C, ^ 4°C above bulk and ~ 16 °C, 8°C, resp., above pore freezing. 
It can be traced to a sudden inversion of the surface tension's T-slope, and thus to a decrease 
in surface entropy at the advancing pore menisci, characteristic of the formation of a single solid 
monolayer of rectified molecules, known as surface freezing from macroscopic, quiescent tetracosane 
melts. The imbibition speeds, that are the squared prefactors of the observed square-root-of-time 
Lucas- Washburn invasion kinetics, indicate a conserved bulk fluidity and capillarity of the nanopore- 
confined liquid, if we assume a flat lying, sticky hydrocarbon backbone monolayer at the silica walls. 

PACS numbers: 47.61.-k, 87.19. rh, 47.55. nb 



When a liquid wets the inner wall of a capillary, it 
forms a concave meniscus and spontaneously invades, a 
phenomenon called capillary filling or spontaneous imbi- 
bition (SI). For channels with diameters of a few nanome- 
ters the interfacial forces driving this process can defy 
gravitational and viscous forces by orders of magnitude 
[1]. Therefore, the study of SI in such extremely spa- 
tially confined geometries has attracted increasing inter- 
est, both from experiment [2 and theory [3^, '4|. Typical 
curvature radii of the advancing menisci are expected to 
be on the order of the pore radius [3 and thus on the 
order of the size of the building blocks of the invading 
liquids. This may alter the equilibrium structural prop- 
erties of the liquids adjacent to the menisci and the flow 
characteristics in their proximity. For example, recent 
molecular dynamics studies show peculiar molecular rear- 
rangements at advancing imbibition fronts of models for 
polymer melts [3 . Moreover, for a pore of a few nm in di- 
ameter any change in the no-slip velocity boundary con- 
dition, by an in- or decreased mobility of the molecules 
in the wall proximity, is expected to alter SI dynamics 
markedly [1|4]. 

Here we present a T-dependent experimental study of 
SI of the linear hydrocarbon n-tetracosane (n-C24H5o, 
denoted C24 in the following) in silica nanopores. C24 
exhibits surface freezing (SF) [5 , that is the formation 
of a single solid monolayer floating on the bulk melt in 
a T-range between the bulk freezing temperature Tf and 
a temperature Tg. Upon SF the C24 molecules are rec- 
tified, parallel-aligned with their long axis along the sur- 
face normal and the center of mass lattice is hexagonal, 
resulting in a 3 nm thick surfactant- like monolayer (see 
illustration in Fig. [T]) . The abrupt onset of SF at Tg al- 
lows one to switch on (and off) the 2D crystalline phase 
by a small T- variation. The goal of our study is to inves- 
tigate whether this peculiar molecular rearrangement is 
detectable by, and how it affects the imbibition dynamics 
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FIG. 1: (color online), (a) Measured normalized imbibition 
speeds vf (definition see text) of C24 in Vycor with 3.4 nm 
(□) and 4.9 nm (o) mean pore radius, resp., in compari- 
son with values calculated under the assumption of meniscus 
freezing (-). Insets: Illustration of a SF-layer at an advancing 
meniscus (T < Tg) and of a sticky, flat lying boundary layer at 
the silica pore wall in the entire T-range investigated, (b) T- 
dependent surface tension of C24 (courtesy of M. Deutsch). 
Insets: Illustrations of SF at a planar C24 surface . The 
dashed lines correspond to extrapolations of the calculated 
imbibition speeds and of the surface tension below Tg in the 
absence of surface freezing, respectively. 



of, C24 in nanopores. Some attention will also be paid 
to possible velocity slippage at the pore walls, and thus 
to deviations from macroscopic SI theory. 

We selected C24 for our study, since it exhibits the 
largest SF T-range of all pure n-alkanes (Tg = 54°C and 
Tf = 50 °C) [5]. As a mesoporous host we have cho- 
sen Vycor glass (Corning glass) from two sample batches 
differing in the mean pore diameter, only. It is a vir- 
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tually pure fused silica glass permeated by a 3D net- 
work of cylindrical, tortuous pores [6]. We cut blocks 
of ~ 10 mm height from the cylindrical rods (cross sec- 
tional area A ^ 37 mm^). After cleaning [7] the samples 
were attached to an electronic balance with a resolution 
of 10~^ mg and were kept in a closed, T-controlled cop- 
per cell with a T- stability of 0.1 °C, see inset of Fig. |2j 
N2 sorption experiments were performed at — 196°C and 
analyzed within a mean-field model for capillary conden- 
sation yielding a mean pore radius of tq = 3.4 nm and 
4.9 nm for samples from the two Vycor batches, resp., 
and an identical volume porosity ^0 =0.31. 

The mass uptake of Vycor as a function of imbibition 
time, m(t), was recorded gravimetrically after bringing 
the Vycor sample in contact with the liquid surface. It is 
plotted, normalized by A, as a function of time t for a se- 
lected set of T's in Fig. [2] We find monotonically increas- 
ing m(t)-curves and the mass uptake rate decreases with 
t. Neglecting the complex pore structure, one can under- 
stand the imbibition dynamics by simple phenomenolog- 
ical considerations: If a liquid wets or partially wets a 
capillary with radius r (contact angle with the pore wall, 
6 < 90°), it is sucked by the capillary pressure, given 
by the Laplace equation, = 2a cos(9)/r. The tube 
geometry calls for a Hagen-Poiseuille type fiow pattern, 
where the viscous drag implies that the volume fiow rate 
is V{t) oc Pc/{r]l{t)) and thus the corresponding mass 
fiow rate is given by rh cx ppc/{r]l{t)), where p is the 
mass density of the liquid and l{t) refers to the height 
of the column that has already been filled up with liquid 
at time t. This length is proportional to m{t). Hence a 
quite simple relation for m{t) and its time derivative is 
arrived at: m{t)m{t) oc p^Pc/r] (x. (? a cos{Q)/r]. It is 
solved by 



m(t) oc p w — cos(O) Vt . (1) 



In a macroscopic SI experiment gravitational forces are 
to be considered, as well. However, for the tiny tubes or 
pores investigated here, the hydrostatic pressure is neg- 
ligible compared to Pc. Consequently the mass uptake 
in a SI experiment should show a Vt-behavior and it is 
sensitive to the fiuid parameters [cr(T), 7^(T), p{T)] and 
the fiuid wall interaction (via the contact angle 6). Re- 
lation ([1]) is known as Lucas- Washburn behavior [9] and 
the square of the second prefactor as imbibition speed, 
Vi{T). Porous media can be considered as complex net- 
works of simple capillaries, which adds an additional t- 
independent geometry factor to relation ([T]). Therefore, 
we expect a V^-behavior for m(t), too. In accordance 
with this statement, the data points of our experiments 
show excellent agreement with V^-fits for the entire SI 
time and temperature range investigated (see Fig. [2| . 
Since we used matrices with identical pore microstruc- 
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FIG. 2: (color online). Specific mass uptake of Vycor for two 
mean pore radii ro due to C24 imbibition as a function of time 
for selected T's close to, but above the bulk freezing temper- 
ature. Solid lines correspond to v^-fits discussed in the text. 
Inset: Raytracing illustration of the thermostated imbibition 
cell employed for the isothermal capillary rise experiments. 



tures for all isothermal SI experiments we can eliminate 
the unknown geometry factor of relation ([T]) by a calcu- 
lation of normalized isothermal imbibition speeds. We 
divide the prefactors necessary to fit the measured m{t) 
curves by the prefactor value at an arbitrarily chosen T 
(here = 60 °C) and arrive at the normalized values 
vf{T), plotted in Fig. [ija) in the vicinity of Tf. They 
are monotonically decreasing with decreasing T, which 
is attributable to the monotonical increase in r]{T) upon 
cooling [10 . Note, however, that at T = = 54°C we 
observe an abrupt change in the T-slope of vf{T). 
In order to analyze this peculiarity we recall the distinct 
signature of SF in the T-behavior of a, known from bulk 
n-alkanes, that is a change from a small negative, above 
Tg, to a large positive T-slope, below Tg [see Fig. [ijb)]. 
It can be understood from elementary surface thermody- 
namics [5|: The T-slope of cr(T) yields information on the 
surface excess entropy: da/dT = —{Ss — S'b), where Ss 
and Sh are the entropies of the surface and bulk, respec- 
tively. The negative slope of bulk C24 for T > Tg is typ- 
ical of an ordinary liquid surface where the molecules on 
the surfaces are less constrained than those in the bulk, 
thus ^g is slightly larger than S'b, yielding d<j/dT < 0. 
SF and its first-order character results in an abrupt re- 
duction of the surface entropy Ss in such a way that Ss 
is smaller than S^ leading to da/dT > 0. Assuming this 
a{T) kink-behavior along with the T-dependency of r] as 
measured with a cone-plate rheometer [10] , we calculated 
normalized imbibition speeds 'L'*^®°(T), plotted as a solid 
line in Fig. [ija). The quantitative agreement with our 
measured changes in the imbibition dynamics is reason- 
able. Note that a change in the n-alkane/silica interac- 
tion at Tg or pore clogging may also be responsible for 
the behavior observed. However, T-dependent contact 
angle measurements on droplets of C24 on planar silica. 



3 



sensitively depending on this interaction, showed 0(T) of 
6±2° (cos 9 ~ 1) and no measurable change upon cross- 
ing Tg. Clogging and the corresponding decrease in v\ 
should be more pronounced in the smaller pores, which 
we also did not observe. Therefore, we feel encouraged 
to solely attribute the distinct change in the SI dynamics 
to a change in cf(T) at the advancing menisci characteris- 
tic of SF. Moreover, the slope difference in cr(T) inferred 
from linear fits of our measured imbibition speeds be- 
low and above Tg of 1.75 ± 0.3 mN/ (m-K) is within error 
margins equal to 1.40 ± 0.1 mN/(m-K) for SF at bulk 
C24 melts. Since this quantity measures the loss of en- 
tropy of surface molecules upon entering the SF state: 
A(da/dT) = dcr/dT(T < Tg) - dcr/dT(T > Tg) = A5, 
our measurements indicate that the surface ordering at 
the advancing imbibition front leads to a loss in entropy 
comparable to the one observed at quiescent bulk surfaces 
and thus, presumably, to a well ordered 2D crystalline 
state. Having in mind the dynamics of the frozen menis- 
cus along the tortuous pore path, the frequent encoun- 
tering of pore junctions with irregular channel geometries 
and the small curvature radii of the menisci, where large 
defect densities of the 2D crystals are expected [JJ^ , this 
may seem surprising. Note, however, the defect forma- 
tion rate is presumably proportional to the monotonically 
in t decreasing meniscus velocity, ~ {).21/ [ns] 1/ns ex- 
pressed here in terms of C24 all-trans length of 3 nm. 
This rate falls below typical crystallization speeds of n- 
alkanes of 0.1 C24 length/ns [12 already after a few 
nanoseconds of elapsed SI time. Therefore, it is plau- 
sibel that the fast crystallization kinetics can heal out 
any occurring defects on t-scales much smaller than de- 
tectable by SI dynamics. 

SF at bulk n-alkane surfaces occurs over a small T-range 
|Tf — Tg| of a few degrees Celsius [5 . However, is 
expected to be significantly shifted downwards in pore- 
confinement [13 . X-ray diffraction experiments on C24 
in 3.5 nm and 4.5 nm Vycor (see [14 ) indicated freez- 
ing of the pore-confined liquid at T = 38°Cand 42 °C, 
only. Thus, our evidence of an unchanged Tg despite 
the sizeable downward-shift of pore freezing suggest that 
pore-confinement may allow one to establish SF over a 
much larger T-range than known of any other free surface 
of bulk n-alkane melts, i.e. Tg — T^ ~ 16 °C. Unfortu- 
nately, we did not find unambigious hints of the Bragg 
peak pattern typical of the SF state in those diffraction 
experiments. The dominating 70% volume fraction of sil- 
ica along with the coincidence of the first maximum of the 
structure factor of both the silica and the nano-confined 
liquid with the position of the dominating SF-Bragg peak 
render its detection in Vycor particularly challenging. 
However, in agreement with our conjecture on the ex- 
tended range of the SF state, this signature has been 
observed in a T-range from Tg down to the pore freezing 
temperature in diffraction studies on silicon nanochan- 
nels partially filled with n-alkanes (10 nm diameter and 



thin silica pore wall corona) p!5^. 

As outlined in Ref. [7 , more detailed insights regarding 
an SI experiment can be achieved by extending Darcy's 
law for the volume flow rate V through a surface area A 
of a porous host. One arrives at the following generalized 
equation for the Lucas- Washburn law: 

".'(.)= -.^ P) 

The tortuosity r characterizes the connectivity and 
meandering of the pores. It could be inferred from ex- 
periments on the self-diffusion of liquids in porous Vycor 
[l6] to r = 3.6 ±0.4, in accordance with simulations of its 
pore morphology JVT . Furthermore, in Eq. ([2| rh refers 
to the hydrodynamic radius over which a parabolic flow 
profile is expected in the pores. It does not necessarily 
have to agree with the mean pore radius ro because of 
either strongly adsorbed, immobile boundary layers, or 
due to velocity slippage at the pore walls O [18]. Us- 
ing the parameters of the pore morphology determined 
from our sorption isotherm measurements (ro = 3.4 nm, 
= 0.31) we overestimate the absolute imbibition speed 
compared to the measured ones for all T's by ^ 30 %, well 
beyond our error margin of ~ 10 % in the SI experiment. 
However, if we assume a hydrodynamic radius rh which is 
reduced by the thickness of approximately 1 hydrocarbon 
backbone, that is 0.4 nm, along with a cr(T) change char- 
acteristic of SF, we arrive at a good agreement between 
measured and predicted m{t) curves for all T's investi- 
gated (see Fig. [2|. More importantly, also for SI in a 
sample of the second Vycor batch with ro = 4.9 nm, the 
assumption of an identically reduced effective hydraulic 
diameter rh = ro — 0.4 nm= 4.5 nm can account for the 
measured SI dynamics, see Fig. |2j The assumption of 
a flat lying, sticky monolayer is corroborated by forced 
imbibition experiments on n-alkanes in Vycor [19] and 
by studies regarding the thinning of n-alkane films in the 
surface force apparatus [20|. Moreover, X-ray reflectivity 
studies indicate one strongly adsorbed, flat lying mono- 
layer of hydrocarbons on silica [21 . The evidence of SF 
in nanopores, presented above, are also supported by a 
study on capillary condensation of hydrocarbons in slit 
geometry [J^. Also in semi-confined, planar, thin film 
geometries on silica surfaces indications of SF have been 
discussed with regard to the complex wetting and freez- 
ing behavior of n-alkanes |M| • 

We reported a first evidence of SF occurring at advancing 
menisci in silica nanopores. This observation testifies a 
remarkable robustness of this archetypical surface order- 
ing transition both upon nanoscale spatial confinement 
and upon self-propelled movement of the interface in a 
rather complex, tortuous pore network. The influence of 
SF on the SI flow dynamics can entirely be accounted 
for by the change in surface tension typical of SF and 
does not additionally slow-down the meniscus movement. 
This finding supports an impressive flexibility of the SF 
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layer similarly as it has been inferred from capillary wave 
spectroscopy The SF layer can be considered as a 

peculiar case of a surfactant. Therefore, our study ex- 
amplifies also how the formation of a surfactant layer at 
nanoscopic menisci can markedly affect the imbibition 
dynamics in nanocapillaries, similarly as it has been ex- 
plored for macropores For the future, we envision 
experiments on the huge variety of systems exhibiting 
SF, which ranges from alcohols, semi-fiuorinated alka- 
nes, and diols to polymers comprising alkyl chains in the 
backbone or as side chains to liquid alloys, liquid crys- 
tals and tetrahedral liquids [27]. The SI dynamics of all 
these technologically important systems may be affected 
by SF. Our study also suggests that nanopore confine- 
ment may allow one to establish the SF state over much 
larger T-ranges than possible at bulk surfaces. There- 
fore, our findings may motivate studies how SF infiuences 
pore crystallization, given the proven impact of SF on 
n-alkane crystallization both in the bulk and in the spa- 
tially confined state [28 . Finally, we believe molecular 
dynamics simulations in combination with synchrotron 
x-ray diffraction experiments and optical birefringence 
measurements on aligned tubular silica channels [29 may 
allow one to elucidate how the SI velocity profiles at the 
advancing menisci will be altered by a SF layer, how the 
SF layer meets the confining silica walls and how the 
available theories and simulations of SF [30] can be ex- 
tended towards non-equilibrium conditions typical of SI 
investigated here. 
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